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ABSTRACT  High-latitude  electric  potential  and  ionospheric  conductance  patterns  are  presented  and  discussed 
for  the  Northern  Hemisphere  during  the  SUNDIAL  period  of  September  23-26.  1986  using  the  Assimilative 
Mapping  of  Ionospheric  Electrodynamics  (AMIE)  technique  of  Richmond  and  Kamide  (1988).  Data  sources  used 
for  the  model  are  satellite  and  ground  magnetometers,  electron  precipitation  instruments,  incoherent  scatter 
radars,  and  ionospheric  coherent  radars  The  period  was  characterized  by  many  substorms.  and  a  wide  variety  of 
instantaneous  patterns  of  electrodynamic  parameters  on  a  hemispheric  scale  are  derived  throughout  this  period, 
of  which  only  a  few  examples  are  displayed  The  entire  set  of  electric  potentials  and  conductances  are  being  made 
available  through  the  NCAR  CEDAR  Data  Base  for  further  analysis  and  utilization  in  simulation  models. 
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ABSTRACT.  High-latitude  electric  potential  and  ionospheric  conductance  patterns  are  presented  and  discussed 
for  the  Northern  Hemisphere  during  the  SUNDIAL  period  of  September  23-26.  1986  using  the  Assimilative 
Mapping  of  Ionospheric  Electrodynamics  (AMIE)  technique  of  Richmond  and  Kamide  (1988).  Data  sources  used 
for  the  model  are  satellite  and  ground  magnetometers,  electron  precipitation  instruments,  incoherent  scatter 
radars,  and  ionospheric  coherent  radars.  The  period  was  characterized  by  mans  substorms,  and  a  wide  variety  of 
instantaneous  patterns  of  electrodynamic  parameters  on  a  hemispheric  scale  are  derived  throughout  this  period, 
of  which  only  a  few  examples  are  displayed  The  entire  set  of  electric  potentials  and  conductances  are  being  made 
available  through  the  NCAR  CEDAR  Data  Base  for  further  analysis  and  utilization  in  simulation  models. 
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1.  INTRODUCTION 

The  high-latitude  distribution  of  ionospheric  electric 
potential  is  a  sensitive  indicator  of  the  coupling 
between  the  solar  wind,  the  interplanetary'  magnetic 
field  (IMF)  and  the  magnetosphere/ionosphere.  It  is 
also  an  important  input  to  simulation  models  of 
ionospheric  and  thermospheric  dynamics  (e.g.  Crow¬ 
ley  et  ai,  1989a,  b ;  Rasmussen  et  a!.,  1986).  We  have 
inferred  the  patterns  of  convection  and  ionospheric 
conductance  every  10  min  during  the  joint  SUN¬ 
DIAL/GISMOS  period  of  September  23-26,  1986 
using  measurements  from  ground  and  satellite-based 
magnetometers,  coherent  and  incoherent  scatter 
radars,  and  satellite  electron  electrostatic  analyzers  in 
the  Assimilative  Mapping  of  Ionospheric  Elec¬ 
trodynamics  (AMIE)  procedure  of  Richmond  and 
Kamide  (1988).  This  method  yields  the  distribution  of 
the  electrostatic  potential  over  the  northern  polar 
region.  The  computed  potentials  are  derived  self- 
consistently  from  different  measurements  of  diverse 


electrodynamic  parameters  and  provide  the  best  esti¬ 
mates  of  the  large-scale  instantaneous  patterns. 

The  AMIE  technique  is  described  in  detail  in 
Richmond  and  Kamide  (1988).  and  has  been  pre¬ 
viously  applied  to  three  other  campaign  periods. 
These  are :  the  GISMOS  (Global  Ionospheric  Simul¬ 
taneous  Measurements  Of  Substorms)  period  of 
January  18-19,  1984  (Richmond  et  ai,  1988.  1990: 
Knipp,  1989),  the  ETS  (Equinox  Transition  Study) 
period  of  September  19.  1984  (Knipp  et  at.,  1989; 
Knipp.  1989).  and  the  period  of  July  23-24,  1983 
(Knipp.  1989).  These  studies  showed  that  realistic 
electric  potential  patterns  could  be  derived.  In  regions 
where  the  inferred  electric  fields  depend  primarily  on 
ground  magnetometer  data,  the  largest  uncertainties 
were  in  the  specification  of  the  ionospheric  conduc¬ 
tance.  These  uncertainties  can  be  reduced  by  the 
inclusion  of  spacecraft  electron  precipitation  measure¬ 
ments  (Knipp  et  al.,  1989 ;  Knipp.  1989). 

The  purpose  of  the  present  paper  is  to  describe 
concisely  the  types  of  data  utilized  to  obtain  the 

Approved  for  public  release 
distribution  unlimi  t.ed. 


Annate*  Geophxucae.  0992-7689.  9006.  J99  10  S  3  00  C  FGS-Giulhitr-Villiirj  _  _ 

_ _  oi  5 _ a 9  n  9  ft 


convection  patterns  for  this  SUNDIAL  period  of 
study,  to  present  a  few  examples  of  patterns  displaying 
noteworthy  characteristics,  and  to  summarize  the 
general  behaviour  of  high-latitude  electrodynamic 
activity  for  the  days  of  September  23-26,  1986. 
Because  of  potential  interest  in  further  utilization  and 
analysis  of  the  time-varying  convection  patterns,  we 
are  making  them  available  now  in  digital  form  to  the 
scientific  community,  even  though  our  own  analysis  of 
them  has  only  begun. 


2.  GEOPHYSICAL  CONDITIONS 

The  SUNDIAL  period  of  September  22-26,  1986  is 
representative  of  solar  minimum  conditions.  The  daily 
10.7  cm  solar  flux  varied  between  68.3  and  69.4 
during  this  period  The  Kp  was  low  (2  +  or  less)  on 
September  22,  and  then  varied  between  6  and  2  on 
September  23-26,  with  the  lowest  values  during  the 
first  part  of  September  25  and  the  latter  part  of 
September  26.  This  variety  of  conditions  is  reinforced 
in  Figure  la.  which  shows  hourly  averages  of  the 
interplanetary  magnetic  field  north-south  component 
(IMF  B.)  for  this  period.  Lower  values  of  Kp  are 
generally  associated  with  positive  IMF  B..  The  quiet 
day  of  September  22.  1986  was  used  as  the  baseline 
from  which  to  measure  the  ground  magnetic  variations 
on  September  23-26.  1986.  Figure  Id  show's  hourly 
averages  of  the  IMF  Bv.  September  23-26.  1986,  was 
generally  a  period  of  IMF  B ,  negative  (and  Bx 
positive).  The  AE  index  was  computed  using  extreme 
variations  in  the  magnetic  north  component  from  the 
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Figure  2 

The  AL  and  A  l'  calculated  every  10  min  from  iht  19  ground 
magnetometer  stations  between  55  and  76‘  for  the  AMIE  stud\  for 
September  25-26.  1986 


19  ground  magnetometer  stations  used  in  the  study 
which  were  located  between  55  and  76  degrees  mag¬ 
netic  north.  The  resulting  AL  and  AU  indices  com¬ 
puted  every  10  min  are  shown  in  Figure  2.  The 
AL  index  in  particular  is  very  active  in  this  period, 
indicating  the  presence  of  many  substorms. 
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Figure  1 

Hourly  IMF  (•)  B.  and (bl  B,  values 
on  September  22-26.  1986  and  the 
times  of  data  coverage  for  the  various 
data  sources  used  in  the  AMIE  run 
for  September  23-26.  1986. 
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3.  STATISTICAL  INPUTS 


The  AM  IE  technique  starts  from  statistical  models  of 
auroral  conductances  and  electric  potential,  and  mod¬ 
ifies  these  in  accordance  with  the  actual  measurements 
for  the  time  of  interest  The  statistical  models  not  only 
provide  a  hemispheric  distribution,  but  also  provide  a 
reasonable  reference  in  regions  where  no  observations 
are  available. 

For  the  present  study,  we  use  the  statistical  auroral 
conductance  model  of  Fuller-Rowell  and  Evans 
(1987).  which  is  parameterized  by  the  10-level  hemis¬ 
pheric  power  index  defined  in  Foster  et  at.  (1986A). 
This  index  is  normally  derived  from  electron  precipi¬ 
tation  measurements  on  board  the  NOAA  series  of 
satellites,  but  instrumental  problems  made  these  data 
unusable  for  most  of  our  SUNDIAL  period  (Evans, 
private  communication).  Consequently,  we  estimated 
the  appropriate  power  level  from  an  empirical  re¬ 
lationship  with  hourly  values  of  the  hemispheric 
power  and  the  AE  index,  where 

P ( GW )  =  4.28  +  0.0624  •  AE  (nT)  . 

The  correlation  coefficient  for  this  relationship  is 
0.71.  We  used  this  formula  with  instantaneous 
AE  values  where  we  subtracted  25  nT  to  account  for 
the  larger  number  of  stations  used  in  our  index 
compared  to  the  regular  index  (Kroehl.  1989).  This 
relationship  is  analogous  to  the  relation  of  power  to 
Kp  published  by  Maeda  et  al.  (1989). 

Observations  are  used  to  modify  the  statistical  auroral 
conductance  patterns,  which  are  combined  with  an 
empirical  model  of  conductance  produced  by  solar 
ultraviolet  radiation  to  yield  the  total  conductance. 
This  model  of  solar  ultraviolet  produced  conductances 
is  a  new  one  that  includes  the  10.7  cm  solar  flux 
variations  as  well  as  season,  time  of  day.  and  magnetic 
latitude.  This  model  is  based  on  Chatanika  radar 
observations  of  electron  densities  (Kroehl  et  al., 
manuscript  in  preparation).  Since  the  September  1986 
period  of  interest  was  near  the  minimum  of  the  solar 
cycle,  it  is  important  to  have  the  conductance  model 
reflect  solar  minimum  conditions,  especially  in  the 
polar  cap. 

Most  previous  studies  using  the  AMIE  technique 
utilized  the  statistical  electric  potential  model  of 
Foster  et  at.  (1986/>),  which  is  also  based  on  the  10- 
level  hemispheric  power  index,  and  was  derived  from 
seven  years  of  Millstone  Hill  radar  measurements.  In 
a  later  study,  Foster  (1987)  combined  S0ndre 
Str0mf]ord  and  Millstone  Hill  radar  ion  drifts  to 
produce  global  electric  potential  patterns  as  a  function 
of  IMF  B,  and  B.  For  the  present  study,  the  initial 
electric  potential  patterns  in  AMIE  use  the  four 
statistical  patterns  of  Foster  (1987)  if  both  IMF 
component  magnitudes  interpolated  from  hourly 
values  are  greater  than  0.5  nT.  For  smaller  magni¬ 
tudes,  the  statistical  patterns  based  on  the  hemispheric 
power  level  (Foster  et  al.,  1986/>)  provide  the  initial 
electric  potential.  The  parameterization  of  the  statisti¬ 
cal  electric  potential  by  IMF  B,  and  B.  is  new  to  the 
AMIE  procedure,  and  was  chosen  because  the  polar 
cap  region  is  not  well  specified  by  data  during  this 


period.  When  the  polar  cap  is  in  darkness,  and  the 
conductivities  are  small,  the  ground  magnetometers 
cannot  easily  sense  the  effects  of  the  electric  fields 
present  there.  We  discuss  briefly  in  Section  5  the 
sensitivity  of  the  results  to  the  choice  of  statistical 
models. 


4.  DATA  INPUTS 

The  data  sources  that  were  incorporated  into  the 
electric  potential  pattern  determination  included  both 
ground  and  satellite  observations.  There  were  32 
ground  magnetometer  stations,  which  are  listed  in 
Table  1.  The  other  sources  are  listed  in  Table  2  and 
include  three  incoherent  scatter  radars,  two  iono¬ 
spheric  coherent  radars,  one  satellite  magnetometer 
and  particle  precipitation  instruments  on  two  satel¬ 
lites.  The  time  coverage  of  these  data  sources  are 
shown  in  Figure  1  above  the  plot  of  B..  Data  that 
covered  an  extended  latitude  range  were  axeraged 
over  about  2.5“  of  magnetic  latitude  in  order  to  be 
compatible  with  the  spatial  resolution  employed  here 
in  the  AMIE  procedure. 


Table  I 

Ground  Magnetometer  Stations 

These  stations  were  used  in  the  AMIE  procedure  in  the  Sl'NDlAL 
campaign  period  of  September  2?-26.  1986  September  22.  1986  was 
used  as  the  quiet  da>  from  which  differences  were  found  for  the 
disturbed  days. 

(D)  =  needed  to  be  digitized 

(*)  =  provided  conductance  with  formulas  modified  from  Ahn  et  al. 
(198?). 


Station 

Mag  Lat 

Mag  Lon 

Resolute  Bav.  Canada 

8-1.1 

304  3 

Mould  Bay.  Canada 

80.6 

263.5 

Cambridge  Bay.  Canada 

77.8 

299.7 

New  AaJesund,  Norway  (D) 

76.9 

114.7 

Baker  Lake,  Canada 

76.2 

320.0 

•Cape  Chelyuskin,  USSR  (D) 

71.6 

174.2 

•Fort  Churchil),  Canada 

70.3 

326.0 

•Barrow.  USA 

69.7 

24S.1 

•Yellowknife,  Canada 

69.6 

294  4 

’Karssarssuaq.  Greenland 

68.9 

44.0 

•Dixon  Island,  USSR  (D) 

68.3 

154.7 

•Poste-de-la-Baleine.  Canada 

68.0 

353.8 

•Leirvogur.  Iceland  (D) 

66.8 

69.6 

•Tixie  Bay,  USSR  (D) 

66.8 

195.5 

•College.  USA 

64. S 

259.6 

•Meanook,  Canada 

62  6 

300.9 

•Cape  Wellen.  USSR  (D) 

62.5 

242  6 

•Glenlea,  Canada 

60  4 

326.8 

Sitka,  USA 

59.8 

277.8 

Ottawa,  Canada 

58.5 

356.2 

Saint  Johns,  Canada 

57.6 

29.1 

Point  Tunguska,  USSR  (D) 

56.5 

162  5 

Victoria,  Canada 

54.1 

292  4 

Magadan,  USSR  (D) 

53.9 

217.4 

Borok.  USSR  (D) 

53.6 

114.2 

Sverdlovsk,  USSR  (D) 

52.7 

132.5 

Fredericksburg,  USA 

50.2 

355.6 

VaJcntia,  Ireland  (D) 

50.1 

71.8 

Boulder.  USA 

49.5 

315.7 

Irkutsk,  USSR  (D) 

46.6 

176.4 

Bav  Saint  Louis,  VSA 

41.8 

338  5 

Memambctsu,  Japan 

37.4 

213  7 
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Table  2 

Radars  and  Satellites 

These  instruments  were  used  in  the  AM  IE  procedure  in  the 
SUNDIAL  campaign  period  of  September '2?-26.  1986. 


Instrument 

Rang< 

e  ir 

i  Magnetic  Coordinates 

Seodre  Stramfjord,  Greenland 

68 

to 

80 

N. 

42.  T 

E 

EISCAT 

62 

to 

72 

N, 

104. 

7  E 

Millstone  Hill.  USA 

50 

to 

TO 

N, 

•IS  to  30  E 

Goose  Bav.  Canada 

65 

to 

77 

N, 

22.0 

E 

SABRE 

62 

to 

65 

N, 

91.3 

E 

DMSP  FT  magnetometer 

50 

to 

90 

N, 

Dear 

10:30-22:30  MLT 

DMSP  F6  electron  detector 

50 

to 

90 

N, 

near 

6:00-16:00  MLT 

DMSP  FT  electron  detector 

50 

to 

90 

N, 

Dear 

10:30-22:30  MLT 

4.x)  Ground  magnetometers 

All  three  components  of  the  measured  magnetic 
perturbation  at  the  ground  are  utilized  in  this  study, 
including  the  vertical  component.  Knipp  (1989)  has 
detailed  how  these  perturbations  are  utilized  within 
the  AM  IE  procedure.  The  magnetic  perturbations  are 
used  primarily  to  determine  the  ionospheric  electric 
field  and  current  patterns,  but  also  used  to  modify  the 
auroral  conductance  model,  as  described  by 
Richmond  and  Kamide  (1988)  and  Richmond  et  at. 
(1988)  based  on  a  modified  form  of  empirical  relations 
presented  by  Ahn  et  al.  (1983).  The  instantaneous 
values  of  the  magnetometer  data  were  used  in  the 
study  without  averaging  since  the  variations  were 
relatively  smooth  in  time. 

4.6)  Incoherent  scatter  radars 

All  of  the  incoherent  scatter  radars  which  operate  in 
the  high-latitude  Northern  Hemisphere  were  operat¬ 
ing  in  the  SUNDIAL  1986  period  Figure  1  shows  that 
Spndre  Strpmfjord  and  EISCAT  started  at  1100  and 
1300  UT,  respectively,  on  September  23,  1986,  after 
the  second  excursion  of  B.  to  large  positive  values. 
The  two  radars  continued  to  operate  until  1810  and 
2350  UT,  respectively,  on  September  26,  1986.  The 
Millstone  Hill  radar  operated  from  2040  UT  on  the 
23rd  to  0530  UT  on  the  25th. 

The  S0ndre  Strpmfjord  radar  scanned  in  geomagnetic 
latitude  every  18  min  and  provided  ion  velocities 
perpendicular  to  the  field  which  were  averaged  in  bins 
of  2.4*  of  apex  latitude.  The  EISCAT  radar  was 
operated  in  the  CP-3-E  mode  which  had  a  scan  cycle 
of  30  min.  The  bisector  measurements  were  converted 
to  perpendicular  and  parallel  components  and  aver¬ 
aged  in  2.3*  latitude  bins.  Millstone  Hill  operated  with 
alternate  low-elevation  azimuth  scans  from  the  east  to 
the  north-west,  and  magnetic-meridian  elevation 
scans  over  a  period  of  24  min.  The  former  were  used 
to  obtain  binned  line-of-sight  velocities  approximately 
perpendicular  to  the  geomagnetic  field.  For  all  of  the 
radars,  all  data  within  the  period  of  one  complete  scan 
centered  around  the  time  of  calculation  were  accepted 
by  the  AMIE  procedure,  though  the  errors  assigned 
to  the  data  were  increased  for  the  data  that  were  not 
close  to  the  time  of  calculation.  The  Spndre 
Str0mfjord  and  EISCAT  radars  also  provided  near¬ 


overhead  ionospheric  conductance  measurements 
that  were  used  to  help  determine  the  auroral  conduc¬ 
tance  distributions. 

4.r)  Ionospheric  coherent  radars 

The  SABRE  (Sweden  and  Britain  Radar  Experiment) 
VHF  coherent  scatter  radar  provided  E  region  elec¬ 
tron  drift  estimates  (Nielsen  ei  al.,  1983).  The  data 
were  separated  into  two  latitude  regions  between  62 
and  65*  magnetic  latitude  and  averaged  over  5  min 
intervals  over  the  entire  longitude  viewing  area  before 
being  converted  to  electric  field  values.  Since  the 
velocity  measurements  saturate  at  the  phase  velocity 
of  the  £-region  irregularities,  these  electric  fields  may 
be  lower  limits  on  the  actual  values.  The  data  were 
available  on  a  continuous  basis,  with  small  signal 
return  indicating  small  electric  fields,  typically  less 
than  10  to  15  mV/m.  Data  within  i  5  min  of  the 
calculation  time  were  included,  with  larger  error  bars 
for  data  5  min  away. 

The  Goose  Bay  HF  coherent  scatter  radar  (Greenwald 
et  al.,  1985)  provides  F  region  irregularity  drifts, 
which  are  similar  to  the  F  region  ion  drifts 
(Ruohoniemi  et  al.,  1987).  A  lack  of  return  signal 
does  not  indicate  small  ion  velocities,  but  only  the 
absence  of  strong  irregularities.  There  were  several 
data  dropouts  in  the  time  period  as  shown  in  Figure  1. 
The  data  were  averaged  every  10  min  over  2‘  bins  in 
magnetic  latitude  and  over  the  entire  longitude  view¬ 
ing  area  before  being  converted  to  electric  fields.  The 
scanning  region  covered  between  65  and  77'  magnetic 
latitude  for  this  period,  but  more  commonly  only 
provided  one  or  two  positions  in  magnetic  latitude. 
All  data  within  15  min  of  the  calculation  time  were 
used.  This  is  the  first  time  these  data  have  been 
included  in  the  AMIE  procedure. 

4 .d)  DMSP  F7  magnetometer 

The  DMSP  F7  satellite  contained  a  magnetometer 
that  was  flown  in  approximately  the  pre-noon,  pre¬ 
midnight  region  over  the  pole  (Rich  and  Gussenho- 
ven,  1987).  As  described  by  Knipp  (1989).  the  gra¬ 
dient  of  the  horizontal  magnetic  perturbations  along 
the  satellite  track  was  derived  and  averaged  over  40  s 
intervals  (about  30u  km  in  distance)  and  used  to 
determine  the  large-scale  field-aligned  current  struc¬ 
ture.  The  gradient  rather  than  the  magnetic  pertur¬ 
bation  itself  was  used  in  order  to  minimize  the 
influence  of  residual  undetermined  crustal  geomagnet¬ 
ic  fields  as  well  as  the  effects  of  large-scale  ionospheric 
currents,  both  of  which  have  relatively  small  gradients 
at  satellite  altitudes  in  comparison  with  the  effects  of 
local  field-aligned  currents.  Any  data  within  20  min  of 
the  calculation  time  were  included  in  the  AMIE 
procedure,  but  with  larger  errors  assigned  to  the 
temporally  more  distant  data. 

4.e)  DMSP  F6  and  F7  precipitating  electron  data 

Rich  el  al.  (1987)  describe  how  auroral  conductance 
enhancements  can  be  derived  from  the  DMSP  electro¬ 
static  analyzers  using  the  formulas  of  Robinson  et  al. 
(1987).  We  use  40  s  averages  obtained  from  the  F6 
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and  F7  satellites,  which  are  in  sun -synchronous  orbits 
with  approximate  equatorial  crossings  at  dawn  and 
dusk  (F6)  and  at  10  :  30  LT  and  22  :  30  LT  (F7).  We 
included  Southern  Hemisphere  data  mapped  to  North¬ 
ern  geomagnetic  conjugate  locations  on  the  supposi¬ 
tion  that  conjugacy  of  precipitation  is  valid  on  the 
300  km  spatial  scale  we  employ  (see  Mizera  et  al.. 
1987).  However,  the  Southern  Hemisphere  measure¬ 
ments  are  assigned  a  50  9?  larger  error  to  account  for 
imperfect  conjugacy. 
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Figure  3 

(a)  The  polar  cap  potential  drop  calculated  every  10  min  from  the 
AMIE  study  of  September  23-26,  1986.  Also  plotted  are  the  potential 
drops  derived  using  the  vB sin’  (S/2)  formula  involving  the  IMF 
Br  B,  and  the  solar  wind  speed  v  (Reiff  and  Luhmann,  1986). 

(b)  The  accompaning  Joule  heating  derived  from  the  Pedersen 
conductance  and  the  electric  field,  corrected  for  expected  additional 
magnitudes  Neutral  winds  are  assumed  to  be  zero. 


5.  RESULTS 

Figure  3a  is  a  plot  of  the  polar  cap  potential  difference 
for  September  23-26.  1986  derived  from  the  AMIE 
technique  and  the  iB  sin?(  0/2)  formula  of  Reiff  and 
Luhmann  (1986)  using  hourly  averages  of  IMF  and 
solar  wind  data.  Here,  t‘  is  the  solar  wind  speed  in 
km  s,  B  is  the  GSM  IMF.  maximized  at  9  nT.  and 
0  =  arcos  (B./B).  The  potential  from  the  IMF  for¬ 
mula  is  composed  of  57  %  of  the  potential  from  the 
present  hour,  29  %  of  the  potential  from  the  past 
hour,  and  14  %  of  the  potential  from  the  hour  before 
that  to  allow  the  ionosphere  to  remember  past  con¬ 
ditions  (Reiff,  private  communication).  The  corre¬ 
spondence  between  the  AMIE  potentials  and  the  IMF 
potentials  is  quite  good  for  this  campaign,  with  the 
larger  disagreements  occurring  after  IMF  data  gaps. 
The  polar  cap  potential  drop  is  anti-correlated  with 
the  IMF  B and  shows  rapid  fluctuations  between  20 
and  HOkV,  with  a  maximum  potential  drop  of 
127.9  kV  occurring  at  2330  UT  on  September  25. 
1986.  Figure  3b  is  the  corresponding  plot  of  the  Joule 
heating  derived  from  the  Pedersen  conductance  and 
the  electric  field.  Included  in  the  total  Joule  heating  is 
the  correction  to  the  electric  field  magnitude  associat¬ 
ed  with  uncertainties  in  the  electric  field,  as  explained 
by  Richmond  el  al.  (1990).  A  maximum  value  of 
443.8  GW  is  found  at  1340  UT  on  September  25. 
1986,  which  can  be  compared  to  a  derived  hemispheric 
power  input  of  61.0  GW. 

Figures  4a  through  4 d  show  the  electric  potential  and 
the  electric  and  magnetic  data  on  September  24  at 
0030  UT  and  0100  UT  when  numerous  radar  electric 
field  data  defined  the  night  side  potential  values  All 
five  radars  were  operating  during  these  times  (M  = 
Millstone  Hill,  G  =  Goose  Bay.  S  =  Spndre 
Slrpmfjord,  A  =  SABRE,  and  E  =  EISCAT).  and 
help  produce  a  distinct  change  in  the  character  of  the 
potential  pattern  near  midnight  between  0030  and 
0100  UT.  Such  coverage  can  be  useful  for  local  studies 
of  substorms  or  other  phenomena. 

There  are  indications  that  the  convection  pattern 
undergoes  major  global  changes  at  times,  in  addition 
to  the  regional  changes  shown  in  Figure  4.  Figure  5 
shows  an  event  where  the  two-cell  convection  appears 
to  be  distorted  by  a  general  clockwise  rotation  and  a 
strong  intensification  of  the  positive  potential  cell  at 
1400  UT  on  September  25,  1986.  Figures  2  and  3  show 
this  to  be  just  after  a  time  when  the  AL  index,  the 
polar  cap  potential  drop,  and  the  Joule  heating  have 
among  the  most  extreme  values  for  the  entire  period. 
In  particular,  the  AL  index  has  been  less  than 
-  100  nT  for  three  hours,  and  is  now  -  860  nT.  This 
value  comes  from  a  station  in  the  midnight  sector, 
indicating  the  presence  of  strong  westward  currents 
associated  with  the  substorm  expansion  current  system 
(Kroehl  and  Kamide,  1985).  The  1  min  IMF  B ,  values 
are  generally  negative  around  this  time,  with  only  2 
brief  excursions  to  positive  values  around  1320  and 
1345  UT.  The  1  min  B.  values  have  been  negative  for 
just  over  1  h. 

In  Figure  5 a,  the  negative  cell  in  the  afternoon  is  fairly 
well  defined  by  the  radar  measurements  for  this  case. 
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Figure  4 

(a)  Electric  potential  derived  for  003U  UT  on  September  24,  1986.  with  superposed  electric  field  observations,  which  have  been  roiohd  90 
counterclockwise  to  show  the  direction  of  plasma  flow .  Contours  are  spaced  at  lOkV.  and  are  dashed  in  regions  where  the  AMIE  prmedurc 
infers  an  uncertainty  in  the  large  scale  electric  field  greater  than  50  9?.  not  including  the  effects  of  uncertainties  in  the  conductance  Tin  letters  M. 
G.  S.  A.  and  E  indicate  the  observations  obtained  by  the  Millstone  Hill.  Goose  Bay.  Sfhidrc  Sirpniford.  SABRE  and  EISCAT  radars, 
respectively.  The  Millstone  Hill  vectors  are  dashed  to  indicate  that  only  the  line-of-sight  velocity  component  M as  used  as  input  data  with  the 
transverse  component  being  derived  from  the  AMIE  fit.  I’eclors  in  the  long-short  dash  pattern  have  large  error  bars  compared  to  other  u<  tore  At 
the  upper  left  are  given  the  total  potential  drop  A <P.  the  initial  hemispheric  power  level  I.  and  the  IMF  B.  and  B:  values  in  nanotcslo  interpolated 
from  hourly  values  (b)  The  ground  magnetic  perturbations  used  to  help  determine  the  electric  potential  shown  in  (a).  Magneto  stations  I, waled 
below  50‘  are  not  plotted,  but  are  included  in  the  calculations.  The  horizontal  perturbations  are  rotated  90'  clockw  ise  to  indicate  the  direction  of 
overhead  equivalent  current,  while  the  vertical  perturbations  are  shown  by  stars  (downward)  or  hexagons  /upward),  (c)  As  in  (a),  but  for 
0)00  UT.  (d)  As  in  (b),  but  for  0100  UT.  and  with  the  addition  of  DMSP  F7  magnetic  perturbation  gradients  along  the  satellite  iratk.  w  hich 
passes  from  the  lop  tow  ard  the  bottom.  Perturbations  lying  on  the  daw  n  side  of  the  track  arc  associated  w  ith  upw  ard  field-aligned  currents  w  hile 
perturbations  tow  ards  the  dusk  side  are  associated  with  downw  ard  field-aligned  currents.  A  scale  of  I  pT,m  is  given  at  the  low  er  le  ft  lor  these 
perturbations 


whereas  the  positive  cell  is  determined  primarily  by 
the  ground  and  satellite  magnetometer  data.  Since  the 
inference  of  electric  fields  from  the  magnetometer 
data  is  sensitive  to  the  conductance  model,  which  has 
its  own  uncertainties,  the  characteristics  of  the  positive 
potential  cel!  in  Figure  5 a  are  less  reliable  than  those 
of  the  negative  cell.  Fortunately,  the  coverage  of 
conductance  data  for  this  particular  time  was  relatively 
extensive  owing  to  the  DMSP  F6  and  F7  overflights. 
Figure  5 d  shows  the  auroral  Hall  conductances  infer¬ 
red  from  these  satellite  data,  along  with  total  (auroral 
and  solar  ultraviolet)  conductances  measured  at  Spnd- 


re  Strpmfjord  and  inferred  from  ground  magnet¬ 
ometers.  superposed  on  the  fitted  distribution  of  Hall 
conductances  determined  by  AMIE.  For  comparison 
the  unmodified  statistical  conductance  model  for  this 
time  is  shown  in  Figure  5c.  As  is  typically  found,  the 
observations  tend  to  make  the  region  of  enhanced 
auroral  conductances  in  Figure  5 d  narrower,  with 
greater  peak  values,  than  given  by  the  statistical 
model  in  Figure  5c.  However,  the  hemispheric  power 
deduced  from  the  fitted  conductance  and  associated 
auroral  energy  flux  patterns  is  somewhat  reduced, 
leading  to  a  lower  power  level. 
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Figure  5 

(a)  4s  in  Figure  4 a.  but  for  1400  UT  on  September  25.  1986  (b)  .4*  in  Figure  4 d.  but  for  1400  UT  on  September  25.  1986  Sole  the  change  in 
scab  of  tin  ground  magnetic  perturbations,  (c)  Statistical  model  of  Hall  conductance  obtained  by  combining  the  auroral  model  of  Fuller- Rowell 
and  Ei ant  (1987 )  for  a  hemispheric  poster  index  of  9  with  a  solar  ultraviolet  conductance  model  appropriate  to  1400  UT  on  September  25.  1986 
Contour  intervals  are  3  Siemens  (d)  Hall  modified  conductance,  with  superposed  data,  for  1400  UT  on  September  25.  1986  Contour  mien  a/s 
are  3  Siemens  The  observations  are  indicated  by  a  vector  extending  towards  the  west  of  the  observing  location,  with  a  scale  shown  at  the  lower 
right  of  the  figure  Thirteen  ground  magnetometer  stations  contribute  total  ( auroral  combined  with  solar  ultraviolet  /  conductance  measurements  in 
the  auroral  zone  The  letters  «  S  »  and  «  E  »  show  the  locations  of  the  S/htdrc  Strpmfjord  and  E1SCAT  measurements  of  the  total  conductance. 
DMSP  F7  auroral  conductances  inferred  from  Northern  Hemisphere  electron  precipitation  measurements  arc  shown  by  solid  vectors  where  a 
•<  o  »  indicates  a  satellite  lime  before  the  calculation  and  a  «  +  »  indicates  a  satellite  lime  afterwards.  DMSP  F6  auroral  conductances  infe  rred  from 
Southern  Hemisphere  data  under  the  assumptions  of  conjugaey  are  shown  by  dashed  vectors.  The  estimated  hemispheric  poster  index  at  this  time 
was  8,  reduced  from  the  initial  level  of  9  by  the  measurements. 


Figures  5e  and  5/ are  plots  of  the  total  horizontal  and 
field-aligned  currents  for  this  time.  The  eastward 
auroral  electrojel  in  the  afternoon  is  associated  with 
the  negative  electrostatic  potential  cell.  The  westward 
electrojet  in  the  morning  and  premidnight  sectors  is 
associated  with  the  positive  cell  which  is  strongly 
distorted  from  a  symmetric  convection-driven  distri¬ 
bution  on  the  nightside  by  auroral  substorm  processes 
(Kroehl  and  Kamide,  1985).  The  field-aligned  currents 
in  Figure  5/ display  well-defined  Region  /  and  Region 
//current  systems  (lijima  and  Potemra.  1978),  where 
Region  /  is  the  downward  current  near  0600  LT  and 
the  upward  current  near  1800  LT.  The  gradient  in  the 
magnetic  perturbation  field  measured  by  the  DMSP 


F7  magnetometer  is  superposed  on  the  field-aligned 
current  distribution.  Since  the  satellite  was  traveling 
from  the  dayside  to  the  nightside.  the  along-track 
magnetic  perturbation  gradients  tend  to  be  associated 
with  downward  field-aligned  currents  when  the  plotted 
vector  is  directed  toward  the  dusk  side  of  the  track, 
and  upward  currents  when  the  vector  is  toward  the 
dawn  side.  Upward  currents  are  usually  associated 
with  the  more  intense  areas  of  particle  precipitation. 
The  extension  of  the  westward  electrojet  and  the 
associated  downward  field-aligned  currents  past  the 
dusk  meridian  is  also  found  in  the  study  of  Kamide  cl 
at.  (1982).  and  appears  to  be  related  to  a  later  phase  in 
substorm  processes.  The  distorted  morning  cell  and 


Figure  Se.  f 

(e)  Height-integrated  horizontal  ionospheric  current  for  1400  UT  on 
Septemher  25.  1986.  w  here  a  vector  of  length  equal  to  10'  in  latitude 
corresponds  to  2  A  m  (0  Field-aligned  current  density  derived  for 
1400  UT  on  September  25,  1986.  Solid  contours  represented  down¬ 
ward  current \  The  contour  interval  is  0.5  jiA  nr  starling  at 
*  0.25  (iA  m'.  The  DMSP  F7  magnetometer  perturbations,  as 
described  in  Figure  Ad,  are  superposed. 


associated  westward  substorm  expansion  current  in 
the  midnight  sector  recurs  several  times  in  the  four 
day  period,  usually  near  or  just  following  peaks  in  the 
AL  (or  AE)  index. 

In  general,  the  convection  for  this  4-day  period  often 
tended  to  have  features  commonly  associated  with 
convection  under  negative  IMF  By  conditions, 
namely,  an  expanded  dawn  convection  cell,  a  stronger 
convection  reversal  shear  in  the  dusk  cell,  and  a 
northeastward  flow  in  the  dayside  throat  region  (de  la 
Beaujardiere  et  at.,  1986;  Foster  el  at.,  1986a; 
Heppner  and  Maynard,  1987).  Figure  5 a  is  one 
example  of  this,  although  it  also  displays  the 
additional  substorm  features  on  the  nightside  men¬ 
tioned  above.  During  periods  when  the  IMF  By  was 
positive,  the  dusk  cell  usually  expanded  relative  to  the 


dawn  cell,  again  in  general  accord  with  previous 
statistical  analyses. 

We  carried  out  several  sensitivity  studies  to  test  the 
stability  and  robustness  of  the  derived  patterns  of 
electric  potential,  current,  and  Joule  heating  to 
changes  in  the  initial  statistical  potential  model 
employed  and  to  deletions  of  various  subsets  of  the 
data.  As  expected  on  the  basis  of  previous  tests 
(Richmond  et  at.,  1988  ;  Knipp.  1989).  noticeable 
changes  in  the  patterns  follow  from  changes  in  the 
inputs  to  the  AM  IE  procedure,  but  most  of  the  global 
characteristics  tend  to  be  preserved  among  the  various 
tests.  The  use  of  an  IMF-parameterized  statistical 
potential  (Foster,  1987)  in  place  of  one  not  based  on 
the  IMF  (Foster  et  at.,  1986 b)  changed  the  shape  of 
the  electric  potential  contours  in  the  polar  cap. 
particularly  when  direct  electric  field  observations  in 
the  polar  cap  were  missing.  We  plan  to  present  in 
more  detail  the  results  of  the  sensitivity  studies  in  a 
future  publication. 


6.  SUMMARY 

The  synthesis  of  several  diverse  data  sets  relating  to 
high-latitude  electric  fields  and  currents  has  enabled 
us  to  obtain  an  approximate  picture  of  the  time- 
varying  convection  patterns  in  the  Northern  polar 
region  for  September  23-26,  1986.  The  ability  to 
include  coherent  scatter  as  well  as  incoherent  scatter 
radar  data  has  improved  the  ability  to  define  the 
convection  in  certain  regions  and  for  longer  periods, 
while  the  inclusion  of  the  three-component  ground 
magnetometer  data  and  horizontal  satellite  magnet¬ 
ometer  data  has  also  helped  to  constrain  the  derived 
patterns.  Other  significant  improvements  we  have 
introduced  include  the  utilization  of  an  IMF  B ,  and 
B.  dependent  statistical  electric  potential  model  as  a 
point  of  departure,  and  a  solar  cycle  dependent  model 
of  the  solar  ultraviolet  component  of  ionospheric 
conductance. 

The  period  September  23-26,  1986  contained  consider¬ 
able  substorm  activity.  During  the  later  phase  of  the 
substorms,  the  dawn  convection  cell  often  expanded 
into  the  pre-midnight  region  in  conjunction  with  the 
westward  electrojet.  IMF  R, -dependent  features  in 
the  relative  sizes  and  orientations  of  the  dawn  and 
dusk  convection  cells  were  often  apparent,  that  gener¬ 
ally  agreed  with  features  deduced  statistically  in 
previous  studies. 

The  electric  potentials  and  conductances  should  be 
useful  for  further  modeling  studies,  and  are  being 
made  available  through  the  NCAR  CEDAR  Data 
Base  in  digital  form  every  10  min  from  September  23 
at  0000  UT  to  September  26  at  2350  UT  on  a  grid  of 
the  Northern  Hemisphere  above  50*  apex  latitude. 
The  grid  spacing  is  every  2*  from  50  to  90  and  every 
hour  of  magnetic  local  time. 
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